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Abstract 

Supercritical branching processes in constant environment conditioned on even- 
tual extinction are known to be subcritical branching processes. The case of random 
environment is more subtle. A supercritical branching diffusion in random environ- 
■ ment (BDRE) conditioned on eventual extinction of the population is not a BDRE. 

However the quenched law of the population size of a supercritical BDRE condi- 



tioned on eventual extinction is equal to the quenched law of the population size of 
a subcritical BDRE. As a consequence, supercritical BDREs have a phase transition 
which is similar to a well-known phase transition of subcritical branching processes 
in random environment. 



1 Introduction and main results 

Branching processes in random environment (BPREs) have attracted considerable interest 
in recent years, see e.g. [21I31IH] and the references therein. On the one hand this is due to the 
more realistic model compared with classical branching processes. On the other hand this 
is due to interesting properties such as a phase transition in the subcritical regime. Let us 
recall this phase transition. In the strongly subcritical regime, the survival probability of a 
BPRE (Z^) t >o scales like its expectation, that is, P(Z^ > 0) ~ const ■ E(Z^) as t — > oo 
where const is some constant in (0,oo). In the weakly subcritical regime, the survival 
probability decreases at a different exponential rate. The intermediate subcritical regime 
is in between the other two cases. Understanding the differences of these three regimes 
is one motivation of the literature cited above. The main observation of this article is a 
similar phase transition in the supercritical regime. 
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Let us introduce the model. We consider a diffusion approximation of BPREs as this is 
mathematically more convenient. The diffusion approximation of BPREs is due to Kurtz 
(1978) and had been conjectured (slightly inaccurately) bei Keiding (1975). We follow 
Boinghoff and Hutzenthaler (2011) and denote this diffusion approximation as branching 
diffusion in random environment (BDRE). For every n G N := {1,2,...} let (Z^ )keNo 
be a branching process in the random environment (Q^, • • • ) which is a sequence of 
independent, identically distributed offspring distributions. If m(Q^) denotes the mean 

offspring number for k G N, then := VnYliIi^°s{ m {Qi n ^))^ k denotes the 

associated random walk where n G N. Set [t\ := max{m G No: m < t] for t > 0. Let 
the environment be such that (S\tn\/V™)t>o conver g es to a Brownian motion (St)t>o with 
infinitesimal drift a G R and infinitesimal standard deviation o~ e G [0, oo) as n — > oo. 
Furthermore assume the mean offspring standard deviation to stabilize at G [0,oo). If 
Zq 1 ^ jn — > z G [0, oo) as n — >■ oo and if a third moment condition holds, then 



in the Skorohod topology (see e.g. |7]) where the limiting diffusion is the unique solution 
of the stochastic differential equations (SDEs) 



for t > where Z = z and S = 0. The processes « J )t>o and (Wn t > are independent 
standard Brownian motions. Throughout the paper the notations F z and ~E Z refer to Z = z 
and So = for z G [0, oo). The diffusion approximation (JTJ is due to Kurtz (1978) (see 
also jS]). Note that the random environment affects the limiting diffusion only through the 
mean branching variance of and through the associated random walk. 

We denote the process (St)t>o as associated Brownian motion. This process plays a 
central role. For example it determines the conditional expectation of Z t 



for every z G [0, oo) and t > 0. Moreover the infinitesimal drift a of the associated 
Brownian motion determines the type of criticality. The BDRE (T2]) is supercritical (i.e. 
positive survival probability) if a > 0, critical if a = and subcritical if a < 0, see Theorem 
4 of Boinghoff and Hutzenthaler (2011). We will refer to a as criticality parameter. 

Afanasyev (1979) was the first to discover different regimes for the survival probability 
of a BPRE in the subcritical regime (see [2TF41IT3] for recent articles). The following char- 
acterisation for the BDRE (j2J) is due to Boinghoff and Hutzenthaler (2011). The survival 
probability of (Z t ) t >o decays like the expectation, that is, ¥(Z t > 0) ~ const ■ E(Z t ) = 
const-exp [(a+^-)tj as t — > oo, if and only if a < —o\ (strongly subcritical regime). In the 




(1) 




(2) 



E z [Z t \(S s ) s < t ] 




(3) 
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intermediate subcritical regime a = —o~ 2 , we have that F(Z t > 0) ~ const ■ t 2 exp {—^ft) 
as t — > 00. Finally the survival probability decays like f(Z t > 0) ~ const ■ t~l exp (—^t) 
as t — > 00 in the weakly subcritical regime a G (—erf, 0). 

This article concentrates on the supercritical regime a > 0. Our main observation is 
that there is a phase transition which is similar to the subcritical regime. Such a phase 
transition has not been reported for BPREs yet. We condition on the event {Z^ = 
0} = {lim^oo Z t = 0} of eventual extinction and propose the following notation. If 
f(Z t > 0|Zoo = 0) ~ const-^{Z t \Z ^ = 0) as t — > 00, then we say that the BDRE (Z t , S t )t>o 
is strongly supercritical. If the probability of survival up to time t > conditioned on 
extinction decays at a different exponential rate as t — > 00, then we refer to (Z t , S t )t>o as 
weakly supercritical. The intermediate regime is referred to as intermediate supercritical 
regime. Our first theorem provides the following characterisation. The BDRE is strongly 
supercritical if a > a 2 , intermediate supercritical if a = a 2 and weakly supercritical if 
a G (0,a 2 ). 

Theorem 1. Assume a,o~ e ,o-b € (0, 00). Let (Z t , S t )t>o be the unique solution of (T5]) with 
Sq = 0. Then 



lim y/fe^ W [Z t > 



t— >oo 



limv^e^'P 2 ^ > 

E— ►£» \ 



lime 

t— >oo 



P 2 > 



Z, 

z, 
z~ 



— — / f(a)4>B(a)da>0 ifae(0,a e ) 



2za e 



-E 



or 



aip(a) da > 
1 



G 



2(^-1) 



> 



if a = o e 
if a > a e 



(4) 
(5) 
(6) 



/or every z G (0, 00) where G v is gamma- distributed with shape parameter v G (0, 00) and 
scale parameter 1, that is, 



F(G U G dx) 



x v e x dx, x G (0, 00) 



(7) 



where ip: (0, 00) — > (0, 00) is defined as 

V2 1 



V(a) := 



7r \ a 



J exp ^ — a( cosh(y)) 2 j cosh(y) c?y, a G (0, 00) 



where f3 :— and where 4>p \ (0, 00) — > (0, 00) is defined as 



</0 



sinh(£) cosh(£)£ 
{u + a{cosh(0) 2 ) {f3+2)/2 



e -a a -m u w-m e -u 



d£ du 



for every a G (0, 00). 



(9) 
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The proof is deferred to Section [2j 

Let us recall the behavior of Feller's branching diffusion, that is, ([2]) with o~ e = 0, which 
is a branching diffusion in a constant environment. The supercritical Feller diffusion con- 
ditioned on eventual extinction agrees in distribution with a subcritical Feller diffusion. 
This is a general property of branching processes in constant environment, see e.g. Jagers 
and Lageras (2008) for the case of general branching processes (Crump-Mode- Jagers pro- 
cesses). Knowing this, Theorem [1] might not be surprising. However, the case of random 
environment is different. It turns out that the supercritical BDRE (Z t , S t )t>o conditioned 
on {Zoo = 0} is a two-dimensional diffusion which is not a BDRE if o~ b > 0. 

Theorem 2. Let cr e £ (0, oo) and cr b , z £ [0, oo) and assume a b + z > 0. If [Z t , S t ) t>0 is 
the solution of ([2]) with criticality parameter a > 0, then 

£ ((Z t , $) t > \Z^ = 0) = £ ((4 S t ) t > ) (10) 
where S t ) t>0 is a two-dimensional diffusion satisfying Z = Z , So = and 

dZ t = (\a 2 e - 2a -S - ) Z t dt + Z t dS t + Jo^¥ t dW t (b) 

a - 2a ° e 1 - ) dt + ^ e dWi e) 
alZ t + a 2 J 

fort> 0. 

The proof is deferred to Section [2j 

It is rather intuitive that the conditioned process is not a subcritical BDRE if > 0. 
The supercritical BDRE has a positive probability of extinction. Thus extinction does not 
require the associated Brownian motion to have negative drift. As long as the BDRE stays 
small, extinction is possible despite the positive drift of the associated Brownian motion. 
Note that if Z t is small for some t > 0, then the drift term of St is close to a. Being doomed 
to extinction, the conditioned process (Z t )t>o is not allowed to grow to infinity. If Z t is large 
for some t > 0, then the drift term of St is close to —a which leads a decrease of (Z t )t>o- 
The situation is rather different in the case o~b = 0. Then the extinction probability of the 
BDRE is zero. So the drift of (St)t>o needs to be negative in order to guarantee Z t — >■ 
as t — )• oo. It turns out that if a b = 0, then the drift of (S t )t>o is —ot and (Z t , S t )t>o is a 
subcritical BDRE with criticality parameter —a. 

We have seen that conditioning a supercritical BDRE on extinction does - in general 
- not result in a subcritical BDRE. However, if we condition (Z t ,S t )t>o on {S^ = — oo}, 
then the conditioned process turns out to be a subcritical BDRE with criticality parameter 
—a. 

Theorem 3. Let a e £ (0, oo) and z £ [0, oo) and assume a b + z > 0. Let {z[ a , St**) t>0 
be the solution of (J2D with criticality parameter a for every a £ R. If a > 0, then 

£ ((zt\S^) t>o l^oo = -oo) = £ Uzt a \st a) ) t> ^ (12) 
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where Zq °^ = Zq*\ 

The proof is deferred to Section [2j 

Now we come to a somewhat surprising observation. We will show that the quenched 
law of (Z t ) t >o conditioned on extinction agrees in law with the quenched law of a subcritical 
BDRE. More formally, inserting the second equation of (fTTj) into the equation for dZ t we 
see that 

dZ t = Qa e 2 - a J Z t dt + a e Z t dW t (e) + yQ^c/Wf } (13) 

for t > 0. This is the SDE for the population growth of a subcritical BDRE with criticality 
parameter —a. As the solution of f JT3|) is unique, this proves the following corollary of 
Theorem [2j 

Corollary 4. Let o~ e G (0, oo) and a^, z G [0, oo) and assume a^ + z > 0. Let {Z\ , 5* ) t>0 
be the solution of (j5J) with criticality parameter a for every a G R. If a > 0, then the 
quenched law conditioned on extinction agrees with the quenched law of the BDRE with 
criticality parameter —a, that is, 

£ ((Z^o^ = 0) = £ [(Zt a) )t> ) (14) 

where Zq~°^ = Zq*\ 

So far we considered the event of extinction. Next we condition the BDRE on the event 
{Zoo > 0} := { lim Z t = oo} of non-extinction. Define U : [0, oo) — > [0, oo) by 

t— s>oo 

U(z):=(a 2 e z + aiy^l (15) 
for z G [0, oo). We agree on the convention that 

oo if c £ (0, oo] c „ . , . 

v J _ := for c G 0, oo) and that ■ oo = 0. (16) 
10 ifc = oo 

Theorem 5. Let a e G (0, oo) and at,, z G [0, oo) and assume a\, + z > 0. Let S t ) t>Q be 
the solution of (J2]) with criticality parameter a > 0. Then 

£ {(Z t ,S t ) t > \Z^ > 0) = £ {(z u S^ t> Jj (17) 
where [Z t , S t ) t>0 is a two-dimensional diffusion satisfying Z = Z , Sq = and 

^ + 2o^ ^1 

2 a*Z t + alU(0)-U(Z t ) 



dZ t = I > 2 + 2a — ^ v ~'' - 1 Z t dt + Z t d£ t + \/ff^ t( |# 



a + 2a ^ . 

a\Z t + al U(0) - U(Z t ) 1 
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(<0 



fort > 0. The quenched law conditioned on non- extinction satisfies that 

£((Z t ) t > |Z oo >0) =c({Z t ) 

where {Zt) t>0 i g the solution of the one- dimensional SDE satisfying Zq = Zq and 



(19) 



dZ, 



-a p + a + 2a 



U(Z t ) 



U(0) - u(z t ) 



Z t dt + a e Z t dWl 



(0 



o-\Z t dW\ h) 



(20) 



fort > 0. 



On the event of non-extinction, the population size Z t of a supercritical BDRE grows 
like its expectation K(Z t \S t ) as t — > oo. 

Theorem 6. Let a e G (0, oo) and <?b, z G [0, oo) and assume + z > 0. Let (Z t ,St)t>o 
be the solution of (j2J) m£/i criticality parameter a G R. T/ien (Z t /e St ^ t>Q is a nonnegative 
martingale. Consequently for every initial value Z = z G [0, oo) there exists a random 
variable Y : Q — > [0, oo) such that 



pSt 



as t — > oo almost surely. 



(21) 



The limiting variable is zero if and only if the BDRE goes to extinction, that is, F Z (Y = 
0) = ¥ z (Z OD = 0). In the supercritical case a > 0, the distribution of the limiting variable 
Y satisfies that 



exp (-AF) 



E 



exp 



(22) 



for all z, A G [0, oo) where G v is gamma- distributed with shape parameter v G (0, oo) and 
scale parameter 1, that is, G u has distribution (J7|). 

The proof is deferred to Section [2j In particular, Theorem [6] implies that Z^ := \imZ t 

t— >oo 

exists almost surely and that Z^ G {0, oo} almost surely. 



2 Proofs 

If <t& = and Zq > 0, then the process (Z t )t>o does not hit in finite time almost surely. 
So the interval (0, oo) is a state space for (Z t )t>o if o"& = 0. The following analysis works 
with the state space [0, oo) for the case 01, > and with the state space (0, oo) for the case 
o~b = 0. To avoid case-by-case analysis we assume cr& > for the rest of this section. One 
can check that our proofs also work in the case = if the state space [0, oo) is replaced 
by (0,oo). 



6 



Inserting the associated Brownian motion (St)t>o into the diffusion equation of (Z t )t>o, 
we see that (Z t ) t > solves the SDE 

dZ t = (a + ^a e 2 ) Z t dt + yfii*Z* dWi e) + dW t {b) (23) 

fortG [0,oo). One-dimensional diffusions are well- understood. In particular the scale func- 
tions are known. For the reason of completeness we derive a scale function for fl23|) in the 
following lemma. The generator of (Z t , S t )t>o is the closure of the mapping Q : Cq([0, 00) x 
R) ->■ C([0, 00) x R) given by 

f a 2 \ d d 1 d 2 

Qf(z, s) := la + -f z—f(z, s) + o^/(*, s) + - (a 2 e z 2 + a 2 z) ^/fo s) 

1 <9 2 <9 2 ^ 

+ 2^V /(Z ' S) + CT ^ /(Z ' S) 
for all z G [0, 00), s6R and every / G Cq([0, 00) x R). 

Lemma 7. Assume a e ,ab,a G (0, 00). Define the functions U: [0, 00) — > (0, 00) and 
V: R -> (0, 00) through 

_2a / 2ft \ 

C/(2) := (a e 2 z + a 2 ) ^ and V(s) := exp ( -— s) (25) 

/or z G [0, 00) and s G R. JTien U is a scale function for (Z t )t>o and V is a scale function 
for (S t )t>o, that is, (U(Z t )) t>Q and (V(St)) t>0 are martingales. Moreover we have that 
gU = and QV = 0. 

Proof. Note that U is twice continuously differentiable. Thus we get that 

■2 



^'• V + *)^S +1 )(*^)" 1 " 1 ^ (26) 
= (a 2 z + al)'^' 1 z (-2a 2 - aa 2 e + Ua 2 + V e 2 j 
= 

for all z G [0, 00). Moreover V is twice continuously differentiable and we obtain that 

/ 2ft \ / 2ft \ -2ft 1 2 / 2ft \ /-2ft \ 2 n 

g exp I s =aexp l s +re e Xp ( s ) =0 (27) 



(2J 



e / \ e 



for all s G R. This shows QU = = QV. Now Ito's formula implies that 



dU(Z t ) = gU{Z t ) dt + u'(Z t ) ■ ^a 2 Z 2 dW t {e) + ^a 2 Z t dW t {b) 

dv(s t ) = gv(s t ) dt + v'(s t ) ^ e dw t {e) 

for all t > 0. This proves that (i7(Z t )) and (^(aS*)) > are martingales. □ 
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Lemma 8. Assume o- e ,o~b,<y G (0, oo). Then the semigroup of the BDRE (Z t ,St)t>o con- 
ditioned on extinction satisfies that 

E^[/(Z,,S,)|Z^0]^ E( " , l™ Z '- S '>' , (29) 
the semigroup of the BDRE (Z t ,S t )t>o conditioned on {S^ = — oo} satisfies that 

^nm,s t )K = -^=^^M (30, 

and the semigroup of the BDRE (Z t , S t )t>o conditioned on {Z^ > 0} satisfies that 

*»> W z„ Sl) |z. > o] = ^Z^-m 12 ^ (31) 

for every z G [0, oo) , s G R, t > and every bounded measureable function f: [0, oo) x R — y 
R. 

Proof. Define the first hitting time T x {jf) := inf{t > 0: r)t = x} of x G R for every 
continuous path 77 G C([0, 00), R)) . As V is a scale function for (St)t>o, the optional 
sampling theorem implies that 

P- (TMS) < 00) = Km r (T_ W (S) < T K( S)) = jto ^gLyS) = v^) < 32 > 
for all s G R and iVeN, see Section 6 in [10] for more details. Thus we get that 
[fiZ^Stp^ = -00] = lim [f(Z t ,S t )\T_ N (S) < 00] 

L 1 iV— >oo 1 

_ Um [/(Z f ,g t )P g '(T_ JV (g)<oo)] 

— (33) 

E^ 5 ) [/(Z,,g f )\/(g f )] 

for all z G [0, 00), s G R and t > 0. The proof of the assertions (|29|) and ( 131]) is analogous. 
Note for the proof of (}3~Tj) that 

P z > 0) = P 2 (lim^ = 00) = lim P z (T N (Z) < T (Z)) = (34) 

for every z G [0, 00). □ 

Proof of Theorm [2J, It suffices to identify the generator Q of the conditioned process. 
This generator is the time derivative of the semigroup of the conditioned process at t — 0. 
Let / G C§ ([0,oo) x R, R) be fixed. Define f z (z,s) := / s (z,s) := 
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fzz(z,s) := jj^f(z,s), f ss (z,s) := jj^f(z,s) and f zs (z,s) := -^f(z,s) for z G [0, oo) and 
s6R. Lemma [8] implies that 

Gf(z,s) 

- lim E(Z ' S) MWiZh, Sh) ~ U(z)f(z, s)} /U(z) _ Q{U ■ f)(z, s) 
/i->o h U(z) 

-2 2 

(® + y) z { u 'f + u f*) (*> s ) + a W (*> s ) + y W«) (*> s ) 



+ ^ K 2 ^ 2 + (U"f + 2C/'L + (*, ^) + ^ (#7. + Uf zs ) (z, s) 



(35) 



(gf( z ,s))U(z) + -^(gU(z))f(z,s) 



for all z G [0, oo) and sGR. Now we exploit that QU = and that ^j(z) = — 2a/(a 2 z + a 2 ) 
for z G [0, oo) to obtain that 

a 2 z 

Gf(z, s) = gf(z, s) - 2azf z (z, s) - 2a^-^ — ^f s (z, s) 

viz + a b 

-a + °i\ zf z (z, s) + 1 (a e V + <x 2 z) s) (36) 

a - 2a A f s (z, s) + \cr 2 J ss {z, s) + a 2 zf zs (z, s) 

aiz + atJ 2 

for all 2 G [0, oo), s G R and all / G Cq([0,oo) x R, R). This is the generator of the 
process (II ip . Therefore the BDRE conditioned on extinction has the same distribution as 
the solution of ([II]). □ 

Proof of Theorm [T], The assertion follows from Corollary H] and from Theorem 4 of 
Boinghoff and Hutzenthaler (2011). □ 

Proof of Theorm [3]. As in the proof of Theorem [2] we identify the generator Q of the 
BDRE conditioned on {Soo = — oo}. Similar arguments as in fl35l) and QV = result in 

Gf(z,s) 

= Sf(z t s) + ^2 (z, s) + a\z (^yf^j (z, s) 

= Qf(z, s) - 2af s (z, s) - 2azf z (z, s) 



-a + J s) - af s (z, s) + — - — —f zz (z, s) + y/„(z, s) + a 2 e zf zs (z, s) 

for all z G [0, oo), s G R and all / G Cq([0,oo) x R, R). This is the generator of the 
BDRE with criticality parameter —a. □ 
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Proof of Theorm [6j Ito's formula implies that 



d- 



e~ St dZ t - e~ St Z t dSt + -e~ St Z t a 2 dt - e~ St Z t a 2 dt 



e~ St ^Z t dt + e~ St Ja 2 Z t dW t [b) + ^e~ St Z t a 2 dt - e~ St Z t a 2 dt 



-(b) 



(37) 



for all t > 0. Therefore [Z t j exp(S t )) t>0 is a nonnegative martingale. The martingale 
convergence theorem implies the existence of a random variable Y : Q — > [0, oo) such that 



^ — »y 

e 5 t 



as t — y oo almost surely. 



(3* 



If a < 0, then = almost surely, which implies y = almost surely. 

It remains to determine the distribution of Y in the supercritical regime a > 0. Fix 
z G [0, oo) and A G [0, oo). Dufresne (1990) (see also [H]) showed that 

f-oo 2 

exp (-as - a e W^ e) ) ds = ^G^. (39) 

Moreover we exploit an explicit formula for the Laplace transform of the BDRE (J2J) con- 
ditioned on the environment, see Theorem 2 of Boinghoff and Hutzenthaler (2011). Thus 
we get that 



E z [exp (-Ay)] = limE 2 

t— ¥OC 



limE 

t— >oo 



exp 



exp 



-A- 



Z, 



limE 2 

t— >oo 



exp ( —A 



s Si 



se[o,t] 



f exp ( - tfjda + exp(-5 t ) 



E 



E 



exp 



exp 



\ / °° exp ( - as - a e Ws e) )ds + \ 

z 



(40) 



^2a/<7§ + J 



This shows 



Letting A oo we conclude that 

2 



p 2 (y = o) = e 



exp 



P Z (^oo = 0) 



The last equality follows from Theorem 4 of [5]. 



(41) 



□ 



Proof of Theorm [5l Analogous to the proof of Theorem [2} we identify the generator Q 
of the BDRE conditioned on {Z^ > 0}. Note that 



-U'(z) 



2a 



U{z) 



U{0) - U{z) a 2 e z + a 2 U{0) - U(z) 



(42) 
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for all z G (0, oo). Similar arguments as in (135]) and QU = result in 

Gf{z,s) 

= Gf(z, s) + (a e V + ^) S ) + ^ z -^M- UXi S ) 

U{z) cr 2 e \ , . , / 2 2a \ 

Q + 2 V(0) + Y J ^ S) + T + ^, + ^(0) -£/(,) J f ^ S) 

for all z G [0, oo), s G R and all / G Cq([0, oo) x R, R). Comparing with (1TB1) . we see that 
^ is the generator of flT8|) which implies f lT7|) . Inserting d^t into the equation of dZ t for 
t G [0, oo) shows that (Z t ) t > solves the SDE ([20]). □ 
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